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Abstract 

The first results on next-to-leading order QCD corrections to production of two Z 
bosons in hadronic collisions in the extra dimension model of Randall and Sundrum 
are presented. Various kinematical distributions are obtained to order a s in QCD by 
taking into account all the parton level subprocesses. We estimate the impact of the 
QCD corrections on various observables and find that they are significant. We also 
show the reduction in factorization scale uncertainty when 0(a s ) effects are included. 
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1 Introduction 



The last missing piece of the standard model (SM), the Higgs boson, remains elu- 
sive to this date, and it is hoped that the Large Hadron Collider ( LHC ) will shed 
light on the mechanism of spontaneous symmetry breaking and discover the Higgs 
bosons. Even if it is discovered there remain fundamental issues, such as the hierar- 
chy problem and others, which make us believe in the existence of some new physics 
beyond the standard model. The LHC which will operate at an enormous centre of 
mass energy (\^S = 14 TeV) promises to be a discovery machine and it is hoped 
that some signals of new physics beyond the SM will be observed. Extra dimension 
models [1^] offer an attractive alternative to the supersymmetry to address the hier- 
archy problem. In this paper we will consider the extra dimension model of Randall 
and Sundrum (RS) [3,4]. There are many important discovery channels at the LHC 
such as 77, ZZ, W + W~, jet production. In the SM the production of two Z bosons is 
suppressed as it begins at the order e 4 in the the electromagnetic coupling and also be- 
cause of the large ZZ production threshold. The two Z bosons can couple to Kaluza 
Klein (KK) gravitons, thus ZZ pairs can be produced through virtual graviton ex- 
change at the leading order. These observations make ZZ production one of the im- 
portant discovery channels. 

At hadron colliders Quantum Chromodynamics (QCD) plays an important role 
as the incoming states in any scattering event are the partons, which are described by 
parton distribution functions (pdf s ). The pdf s depend on the factorization scale (/j, f ) 
which is, to a large extent, arbitrary. This scale fip enters into any observable and 
makes it sensitive to the choice of its value and any leading order computation suffers 
from this sensitivity. However, as a computation beyond the leading order (LO) is 
carried out, the fi F dependence partially cancels yielding results less sensitive to the 
factorization scale. It also improves upon LO results in that it includes missing higher 
orders terms of the perturbation series which can be large. It is, thus, the motivation 
of this paper to consider production of Z boson pairs at the LHC at next-to-leading 
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order (NLO) accuracy in the strong coupling constant in RS model. 

Leading order studies for ZZ production in the SM can be found in [5]. Z pair 
with a large transverse momentum jet at LO was studied in [6]. LO study for ZZ 
production in the context of extra dimension model of RS was carried out in [7] and 
coupling of radion with gluon and top quark loop was cosidered in [8]. Because 
of its importance ZZ production has also been studied to NLO accuracy in the SM 
[9-11]. These results were subsequently updated in [12,13]. NLO studies in SM via 
gluon fusion were carried out in [14,15]. These studies provide the precise estimate 
of higher order effects through K factor as well as the sensitivity of the predictions 
to factorization scale. Importantly, the corrections turned out to be larger than the 
expectations based on soft gluon effects justifying a full-fledged NLO computation 
taking into all the processes. We presented NLO results for ZZ production at the 
LHC in large extra dimension model of Arkani-Hamed, Dimopoulos and Dvali [1,2] 
in [16] where it was shown that the K factors are large. Although NLO results are 
available in SM and ADD model they do not exist in literature in the context of RS 
model which is the material of the present paper. 

The results which are presented in this paper are obtained using our NLO Monte 
Carlo code (which is implemented on FORTRAN 77) that can easily accommodate any 
cuts on the final state particles and obtain various kinematical distributions. Our code 
is based on the method of two cutoff phase space slicing (for a review of the method 
see [17]) to deal with soft and collinear singularities in the real emission contributions. 
We will use the matrix elements presented in [16] and refer the reader to this paper 
for further details. 

In what follows we will first briefly describe the RS model and then present the 
numerical results and finally conclude. 
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2 RS Model 



In the RS model the single extra dimension is compactified on a S 1 /Z 2 orbifold 
with a radius R c which is somewhat larger than the Planck length. Two 3-branes, the 
Planck brane and the TeV brane, are located at the orbifold fixed points = 0, n, 
with the SM fields localized on the TeV brane. The five-dimensional metric, which is 
non-factorizable or warped is of the form 

ds 2 = e' 2!CR ^ri^dx^dx u + R 2 c d<j) 2 (1) 

where < <f) < tt. The the huge ratio ~ 10 15 can be generated by the exponent 
JCR C which needs to be only of O(10) thereby providing a way of avoiding the hier- 
archy problem. It was shown in [18, 19] that R c can be stabilized against quantum 
fluctuations by introducing an extra scalar field in the bulk. 

The tower of massive Kaluza-Klein (KK) excitations of the graviton, hffl, interact 
with the SM particles by: 

i ttICR 00 

An, ~ -=T^(x)h^(x)-^^T^(x)h^(x). ( 2 ) 
Mpi Mpi n=1 

is the symmetric energy-momentum tensor for the SM particles on the 3-brane, 
and Mpi is the reduced Planck scale. The masses of the hffl are given by 

M n = x n Ke^, (3) 

where the x n are the zeros of the Bessel function J 1 (x). The zero-mode couples weakly 
and decouples but the couplings of the massive RS gravitons are enhanced by the ex- 
ponential e nKRc leading to interactions of electroweak strength. Consequently, except 
for the overall warp factor in the RS case, the Feynman rules in the RS model are the 
same as those for the ADD case [20,21] 

The basic parameters of the RS model are 

c = KfWi, (4) 
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where m is a scale of the dimension of mass and sets the scale for the masses of the KK 
excitations, and c is an effective coupling. The interaction of massive KK gravitons 
with the SM fields can be written as 

oo 

4„t~-— ^T>»(x)h$(x). (5) 

n 

Since K, is related to the curvature of the fifth dimension we need to restrict it to 
small enough values to avoid effects of strong curvature. On the other hand K, should 
not be too small compared to M P i because that would reintroduce hierarchy. These 
considerations suggest 0.01 < Co < 0.1. For our analysis we choose to work with the 
RS parameters Co and Mi the first excited mode of the graviton rather then m . 
Let us define T>(Q 2 ) as the sum of KK graviton propagators 

°° 1 A 

V{Q2) = ^Q^M^M^^' (6) 

where M n are the masses of the individual resonances ( see Eq. [3] ) and the T n are the 
corresponding widths. The graviton widths are obtained by calculating their decays 
into final states involving SM particles. A is defined as 

_! - m ° ? (7) 

n=l X * X ™ ^ m Xn 

where x s = Q/m . We have to sum over all the resonances to get the value of X(x s ). 
This is done numerically and for a given value of x s , we retain all resonances which 
contribute with a significance greater than one per mil, and treat the remaining KK 
modes as virtual particles (in which case the sum can be done analytically). 

As the gravitons couple to Z bosons, PP — > ZZ can now also proceed through a 
process where gravitons appear at the propagator level. These new channels make it 
possible to observe deviations from SM predictions if extra dimensions exist. In the 
following we will consider spin-2 gravitons only at the propagator level and investi- 
gate this process at NLO level. 
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3 Numerical results 



In this section we present invariant mass (Q) and rapidity (Y) distribution of the Z 
boson pairs. These kinematical variables are defined as 

Q 2 = (p Zl +Pz 2 ) 2 , F = ^ln^, (8) 

z f 2 ■ q 

where Pi and P 2 are the momenta of colliding hadrons, and q = pz 1 + pz 2 denotes the 
sum of the Z-boson 4-momenta. In obtaining these distributions all order a s contri- 
butions have been taken into account. At leading order in SM, the process proceeds 
through qq initiated process. As the gravitons couple to the Z-bosons, qq and gg ini- 
tiated processes with virtual gravitons also contribute at the same order in QCD. We 
have considered all the qq, gg initiated one loop virtual and, qq,qg,gg initiated real 
emission corrections to these processes, both in the SM and the gravity mediated pro- 
cesses, and their interferences. At the virtual level we used method of Passarino and 
Veltman [22] to reduce tensor loop integrals to scalar integrals. In dealing with real 
emission contributions we have used two cutoff phase space slicing method. Here, 
using two small dimensionless slicing parameters 5 S and 5 C , the singular (soft and 
collinear) regions of phase space are separated from the finite hard noncollinear re- 
gion. We will refer to the sum of contributions to crosssection from virtual, soft and 
collinear regions as 2-body contribution, and from hard noncollinear region as 3-body 
contribution. The soft singularities cancel between real and virtual contributions and 
the collinear singularities were removed by mass factorization in MS scheme, this 
gives the finite 2-body contribution. Finally the kinematical distributions were ob- 
tained by integrating the 2-body, 3-body and leading order pieces over the phase space 
using monte carlo methods. Individually 2-body and 3-body contributions depend on 
the slicing parameters 5 S and 5 C but the sum is invariant against variations of these pa- 
rameters over a wide range. For further analysis we will use 5 S = 10 ~ 3 and 5 C = 10~ 5 . 
For further details please see [16]. 

Below we present various distributions for the LHC with a center of mass energy 

5 



of 14 TeV as a default choice. However we will also present some results for a center 
of mass energy of 10 TeV for the LHC. For numerical evaluation, the following SM 
parameters [23] are used 

M z = 91.1876 GeV, sin 2 6 W = 0.231 (9) 

where 9w is the weak mixing angle. For the electromagnetic coupling constant a 
we use a -1 = 128.89. CTEQ6 [24,25] density sets are used for parton distribution 
functions. 2-loop running for the strong coupling constant has been used, and we 
have done calculation with 5 quark flavors and taken the masses of quarks to be zero. 
The value of Aq C d is chosen as prescribed by the CTEQ6 density sets. At leading order 
we use CTEQ6L1 density set ( which uses the LO running a s ) with the corresponding 
Aqcd = 165 MeV. At NLO we use CTEQ6M density set ( which uses 2-loop running 
a s ) with the Aq C d = 226 MeV; this value of Aq C d enters into the evaluation of the 
2-loop strong coupling. The default choice for the renormalization and factorization 
scale is the identification to the invariant mass of the Z boson pair ie., ji F = ji R = Q. 
Furthermore the Z bosons will be constrained to satisfy \yz\ < 2.5, where yz is the 
rapidity of a final state Z boson . 

In Fig. [T] we have plotted the invariant mass distribution both for the SM and the 
signal for LHC at lATeV. The two curves with peaks correspond to the signal and the 
remaining two curves give SM predictions. Here we have chosen c = 0.01 and Mi = 
1500GeV. To highlight the importance of QCD corrections we have also displayed 
the LO results of SM and the signal, and we observe that the K factors (defined as 
K = da NLO /da LO ) are large. For the signal the K factor is 1.82 at Q = 1500 GeV. Also 
note that the gravitons appear as ever widening peaks, these resonance peaks are clear 
signals of the RS model as opposed to the enhancement of the tail in ADD model. 
Next we present in Fig. |2] the effects of varying the parameter cq on the invariant 
mass distribution. All the curves shown correspond to NLO results, and we have also 
plotted the SM background for comparison. 

In Fig. |3] we have plotted the rapidity distribution da/dY at LO and NLO both 
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for SM and the signal for Co = 0.01. We have plotted this distribution in the interval 
—2.0 < Y < 2.0 and have carried out an integration over the invariant mass interval 
1450 < Q < 1550 to increase the signal over the SM background. As expected the 
distribution is symmetric about Y = 0. In Fig. H] we have shown the variation of 
rapidity distribution with c . 

We have mentioned before that the NLO QCD corrections reduce the sensitivity 
of the cross sections to the factorization scale up; this we now show in the Fig. [5] We 
have plotted SM and the signal both at LO and NLO, and have varied the factorization 
scale lip in the range Q/2 < fi F < 2Q. The central curve in a given band (shown by 
the dotted curves) correspond to fi F = Q. In all these results the renormalization scale 
is fixed at ijlr = Q. We notice that the factorization scale uncertainty at NLO is much 
reduced compared to the LO. For instance at Q = 1500 GeV varying fi F between Q/2 
to 2Q shows a variation of 20.6% at LO for the signal, however the NLO result at 
the same Q value shows a variation of 7.1%. Similarly we show the dependence on 
factorization scale at LO and NLO in the rapidity distribution in Fig. [6] and Fig. for 
SM and signal respectively. 

At the end we present in Fig. [8], da/ dQ for LHC with a centre of mass energy of 
10 TeV at NLO both for SM and signal. For comparison we have also plotted the 
14 TeV results in the same figure. 

4 Conclusions 

In this paper we have carried out a full NLO QCD calculation for the production 
of two Z bosons at the LHC at lATeV in the extra dimension model of Randall and 
Sundrum. Here we take all order a s contributions, both in the SM and in the grav- 
ity mediated processes and their interferences, into account. We have presented in- 
variant mass and rapidity distributions both at LO and NLO. We use CTEQ 6L1 and 
CTEQ 6M parton density sets for LO and NLO observables, respectively. Significant 
enhancements over the LO predictions are observed. The K factors are found to be 
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large, for instance in invariant mass distribution the signal has a K factor of 1.82 at 
Q = 1500 GeV (the position of first resonance). We have also presented the effects of 
variation of parameter c both in Q and Y distributions. We have shown that a sig- 
nificant reduction in LO theoretical uncertainty, arising from the factorization scale, 
is achieved by our NLO computation. For instance at 1500 GeV varying fi F between 
Q/2 to 2Q shows a variation of 20.6% at LO for the signal, however the NLO result at 
the same Q value shows a variation of 7.1%. Thus our NLO results are more precise 
than the LO results and suitable for further studies for constraining the parameters 
of the RS model. Invariant mass distribution is also presented for LHC at a center of 
mass energy of lOTeV at the NLO level. 
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Figure 1: Invariant mass distribution for SM and signal both at LO and NLO. Dash- 
dot curves represent LO results and solid curves give NLO results. We have chosen 
Mi = 1500 GeV and the parameter c = 0.01. 
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c variation, LHC (V S =14 TeV) 




Figure 2: Effect of variation of c on invariant mass distribution. All the curves corre- 
spond to NLO results with Mi fixed at 1500 GeV. The solid curve corresponds to SM 
and the dash-dot curves to the signal. The signal is plotted for c = 0.01, 0.04, 0.08 and 
the dash size increases with increasing c 
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LHC (V S =14 TeV) 




c =0.01, M 1= 1500 GeV Y 



Figure 3: Rapidity distribution for SM and signal both at LO and NLO. Dash curves 
represent LO results and solid curves give NLO results. We have chosen Mi = 
1500 GeV and the parameter Co = 0.01. To enhance the signal we have integrated 
over Q in the range 1450 < Q < 1550. 
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c variation, LHC (V S =14 TeV) 




Figure 4: Rapidity distribution for SM and signal at NLO. Dash curves represent sig- 
nal and solid curves gives SM contribution. We have chosen Mi = 1500 GeV and 
plotted the signal for c = 0.01, 0.04, 0.08. To enhance the signal we have integrated 
over Q in the range 1450 < Q < 1550. 
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p. F variation, LHC (V S =14 TeV) 




Figure 5: Factorization scale variation in the invariant mass distribution. The curves 
correspond to c = 0.01 and M x = 1500GeV at the LHC at y/S = 14 TeV. The fi F is 
varied between Q/2 and 2Q. The dash curves correspond to fip = Q 
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H F variation, LHC (V S =14 TeV) 




Figure 6: Factorization scale variation in the rapidity distribution for the SM. The 
curves correspond to c = 0.01 and M x = 1500GeV at the LHC at y/S = 14 TeV. To 
enhance the signal we have integrated over Q in the range 1450 < Q < 1550. The [if 
is varied between Q/2 and 2Q. The dash curves correspond to ji F = Q 
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H F variation, LHC (V S =14 TeV) 
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Figure 7: Factorization scale variation in the rapidity distribution for signal. The 
curves correspond to c = 0.01 and M 1 = 1500GeV at the LHC at VS = 14 TeV. 
To enhance the signal we have integrated over Q in the range 1450 < Q < 1550. The 
pLp is varied between Q/2 and 2Q. The dash curves correspond to fx F = Q 



15 



LHC (V S =14 TeV, 10 TeV) 
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Figure 8: Invariant mass distribution for SM and signal at \/S = WTeV and 14TeV. 
All the curves correspond to NLO results. We have chosen M 1 = 1500 GeV and the 
parameter c = 0.01. 
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